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The effects of the finite top quark width in the top pair production cross section close to the threshold are 
discussed in this talk. We introduce a cross section with a cut on the invariant masses of the top and antitop that 
can be calculated theoretically with effective field theory (EFT) methods. The matching procedure to implement 
the physical phase-space boundaries in the NRQCD framework ( "phase-space matching" ) is briefiy outlined. 



The measurement of the a{e^e tt) line 

shape at a linear collider (LC) operating at en- 
ergies around the top-antitop threshold ^ 
350 GeV) would allow for very accurate determi- 
nations of the mass, the width and the couplings 
of the top quark [T] . In view of the accuracy ob- 
tainable at the LC the theoretical uncertainties 
for the cross section predictions should be low- 
ered to a level of a few percent [2] . 

Close to threshold the top quark pairs are pro- 
duced with small velocities v <^ 1. In the stan- 
dard QCD perturbative expansion singular terms 
~ as/v, as log V arise from the ratios of the phys- 
ical scales involved: the top mass nit, the relative 
three-momentum p rritv and the nonrelativistic 
energy E ~ rritv^ . The summation of these terms 
can be achieved systematically by means of non- 
relativistic QCD (NRQCD), the low-energy EFT 
that describes the quantum fluctuations of full 
QCD for the kinematic situation of top quarks 
close to threshold. Renormalization group im- 
proved calculations of the tt total cross section 
with NNLL order accuracy for the QCD effects 
have become available in this framework [3]. 

A study of electroweak corrections in tt produc- 
tion at threshold beyond leading order has not 
been systematically undertaken until recently [H 
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[S]. The large top width, being of the same or- 
der as the nonrelativistic energy, is essential in 
the description of the tt threshold dynamics. The 
dominant top decay mode in the Standard Model 
is t ^ bW'^, which gives Ft ^ 1.5 GeV. It was 
shown |6j that in the nonrelativistic limit the top- 
quark width can be consistently implemented by 
the replacement E E + iFf in the results for 
the total cross section for stable top quarks. This 
replacement rule accounts for the LL electroweak 
corrections to the total tt cross section, but a com- 
plete treatment at NLL and NNLL requires the 
use of an extended NRQCD formalism accounting 
for the top quark instability. 

In the NRQCD framework the quark width ef- 
fects are incorporated at leading order through 
the bilinear quark operators 

^^ = E^p^r,v^p-i-5]x^^r,xp. (1) 
p p 

They represent a LL effect as compared to the 
typical energy of a Coulombic top quark, of the 
order of a few GeV. The terms in Eq. ^ describe 
an absorptive on-shell process that has been in- 
tegrated out from the theory. In this way, the 
EFT does not contain any details and differen- 
tial information of the decay, but can predict the 
total cross section, an inclusive observable. We 
can use the optical theorem to compute the total 
cross section from the imaginary part of the e^e~ 
forward scattering amplitude (see Fig [T|). The 



Figure 1. Tree-level double-resonant diagrams in 
the full and effective theories. The dashed line 
means we extract the imaginary part or, equiva- 
lently, that we perform the phase space integra- 
tion along the cut lines. 

coefficients in Eq. (jlj are determined from the 
matching procedure of the EFT to the full the- 
ory (QCD and the electroweak theory). The con- 
tributions from real bW final states are thus in- 
cluded through matching conditions of the EFT. 
The LO bilinear quark terms in the NRQCD La- 
grangian lead to the propagator 



(2) 



Po - pV(2TOt) + irt/2 



for a top/antitop with momentum {po,p) ^ 
{mtv'^,mtv). The inclusion of a finite width in 
the top propagator modifies the high-energy be- 
haviour of the phase-space integrations in the 
EFT. Let us illustrate this with an example. Con- 
sider the double-resonant diagram of Fig. [1] in 
the full theory with center of mass momentum 
q — (2mt + E,0). When the quark lines are close 
to their mass shell the 4-particle phase space in- 
tegration reduces to 



/ + 00 /' + OC 
dpo / d|p|p^ — 
-oo Jo \tl 



Ih + ief \t2 + iei^ 
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where the variables ti^2 = ii,2(pojP) = 2rnt(-| 

2 

Po — ^jj) retain the leading order terms in the 
nonrelativistic expansion of the top and anti- 
top off-shellness. The Fj factor in the numer- 
ator arises from the phase-space integration of 
the and bW^ subsystems because in the 

nonrelativistic limit the invariant masses of the 
top/antitop quark can be set to m^. 

We notice that though the phase-space in the 
full theory is cut-off by the top mass, the nonre- 
lativistic limit formally makes nit — ^ oo and the 
phase-space boundaries are thus taken to infin- 



Figure 2. (e+e^)(e+e^) scattering amplitude 
with a kinetic energy correction and a Coulomb 
potential. Its imaginary part gives a NNLO di- 
vergent contribution to au. 

it}0. The corresponding NRQCD amplitude (di- 
agram on the right in Fig. [1]) reproduces Eq. ^ 
when Ft is treated as an insertion. The NRQCD 
powercounting, however, tells that the Tt ^ E 
term needs to be resummed as part of the top 
propagator (see Eq.(2)), thus effectively replac- 
ing replacing ie iTt, to finally give (x v/m^, 
with V = y/{E + iY^^)/mt. 

Despite the integration limits, the tree- level 
integration in Eq. ([3]) is finite. However the 
integrand becomes more sensitive to large mo- 
mentum regions once we include relativistic cor- 
rections ~ p^/mj. Using dimensional regular- 
ization these subleading contributions can lead 
to 1/e singularities if the high energy behaviour 
of the EFT phase-space integration is logarith- 
mic. This occurs for some of the NNLO ma- 
trix elements containing a single insertion of the 
Coulomb potential and a p'^/mf correction. Fig- 
ure [2] shows an example. Similar UV divergences 
had also been noted in Ref. [4], where the imag- 
inary parts of the NRQCD matching conditions 
at NNLL coming from Wb cuts in the full the- 
ory were determined and inserted into NRQCD 
matrix elements. These UV divergences would 
not exist for a stable top quark. These phase 
space divergences can be handled with the usual 
EFT renormalization techniques. In this spe- 
cific example the UV divergences are absorbed 
by (e+e~)(e+e~) operators that will thus have 
complex Wilson coefficients Cv/a{i^) The 
imaginary parts of these Wilson coefficients con- 
tribute to the tt cross section through the opti- 
cal theorem, AtTft — Im[C'v' + Ca\- The renor- 
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malization procedure of course does not deter- 
mine the matching conditions at the hard scale, 
C'v'/a(/^ — "^t); ^ matching prescription is re- 
quired in order to put back into the EFT the in- 
formation about the physical phase-space bound- 
aries. A possible approach is to define the match- 
ing condition as the difference between the full 
theory computation for the hW^hW~ total cross 
section (considering h and W as stable) and the 
NRQCD result, similarly to the procedure used in 
Ref. [7] for the case of W'^W^ threshold produc- 
tion. Obviously, the full theory computation with 
the 4-particle final state is a rather complicated 
task beyond tree-level, even if an asymptotic ex- 
pansion technique is employed. 

At this point, a different path can be followed, 
that will reduce the amount of full theory input 
required and allow us to make contact with the 
experimentally measured cross section. Among 
the experimental cuts that need to be applied to 
select bW^bW" events in tt production, cuts on 
the top and antitop invariant masses will be es- 
sential. Close to the threshold, it is more natural 
to define the cut on the the top invariant mass 
squared pj with respect to its on-shell value, i.e. 

b?-m?|<A^ (4) 

and similarly for the antitop. In the nonrelati- 
vistic limit, the top and antitop off-shellness are 
characterized by the variables ii,2(PoiP)- The 
cutoff A divides the phase-space into three dif- 
ferent regions: the double-resonant region, with 
both |ii^2| < A^, the single-resonant regions, cor- 
responding to (|ii| < A^, 1^2! > A^) or (|t2| < A^, 
1^1 1 > A^), and the non-resonant region, where 
both \ti,2\ > A^. The regions are shown in the 
phase-space diagram in Fig. [3] (gray, light-gray 
and white, respectively). The red lines in Fig. [3] 
correspond to the top and antitop on-shell con- 
ditions ti_2 = 0. For the case Ft = 0, the cut 
through the top and antitop EFT propagators 
yields a delta function instead of a Breit-Wigner 
distribution, and the phase-space shrinks to the 
point given by the intersection of the two on-shell 
lines (IpI = \fmE). Restricting the phase-space 
boundaries to the double-resonant region defines 
a new threshold tt cross section, CTtt(A), that 
corresponds to taking into account only those 
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Figure 3. Phase-space integration region spanned 
in the (po,p) variables. We have chosen E = 
5 GeV and A = 80 GeV for this example. 

measured bW'^ and bW~ events with invariant 
masses inside the window given by the condition 
([!]). The cutoff captures the double-resonant re- 
gion where the top and antitop are close to their 
mass-shell and the NRQCD expansion applies. It 
excludes those parts of the phase space that are 
not described well by the EFT. The latter is true 
provided that the value of A formally satisfies 
(£',r) ^ A^/mt ^ mj. The physical meaning 
of the new matching condition, 

lmCv/A=^u{^)-^fr . (5) 

becomes thus clear in this approach: it corrects 
for the wrong behaviour of the effective propaga- 
tors outside the nonrelativistic domain (ct^^"^ is 
the standard NRQCD result, i.e. for A — > oo). 
We can go a step further by realizing that in- 
side the double-resonant region the full theory 
cutoff cross section Ctf (A) can be nonrelativisti- 
cally expanded since A <C to*. The computation 
of cTff (A) can thus be carried out using NRQCD 
Feynman rules and, in this way, the whole ap- 
proach is interpreted as a matching procedure 
within the EFT itself. There are of course con- 
tributions with the same bW^bW~ final state 
that are not described by the EFT computation 
(the interferences of amplitudes containing a sin- 
gle top/antitop line, for example). The latter 
are part of the irreducible background, whose nu- 
merical contribution to (Ttf (A) can be estimated 
using the tree-level predictions since there is no 
Coulomb-enhancement for these terms at higher 
orders in Ug. 

The coefficients Cy/A provide the corrections to 
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Figure 4. 0(0;^) phase-space matching. 

the standard NRQCD result due to phase-space 
effects at a practical level. From a formal point 
of view, we have overlooked the fact that the cor- 
rections CvjA defined through Eq. (O are energy 
dependent, and cannot be cast as the Wilson coef- 
ficients of local (e"'"e~)(e+e~) operators. A local 
action is needed in order to use the optical theo- 
rem to determine the cross section from the imag- 
inary part of the e^e~ forward scattering ampli- 
tude. The natural solution is to write the correc- 
tions Cy/^ as the coefficients of an operator prod- 
uct expansion in (e+e~)(e~''e~) operators with in- 
creasing powers of time derivatives id^ ^ E. This 
is achieved practically by expanding the phase- 
space effects in powers of E/A. For example, the 
matching of the 0(a°) diagram shown schemati- 
cally in Fig. [4] yields the phase-space correction: 

As expected, the result vanishes for Ff 0. We 
also note that the non-local contribution cx w in 
the diagrams of Fig. 0] cancels in the difference. 
The term proportional to Ff/A gives the domi- 
nant phase-space correction to the tt cross sec- 
tion, whose relevance will depend on the chosen 
value for the cutoff. We mentioned above that 
our approach is formally valid for A^/mt ^ mt- 
This condition may seem necessary due to the 
appearance of powers of A^/m^ when higher di- 
mensional operators are considered. However, a 
numerical analysis of the coefficients of such con- 
tributions shows that the OPE series has a good 
convergence even for values of A close to nit , al- 
though smaller values (A < 0.6mt) are required 
to lower their size below the desired accuracy of 
1%. We shall use the counting A ~ toj to iden- 
tify the order at which the phase-space matching 
corrections contribute. In this scheme, the cor- 
rection from Eq. ^ scales as F/A ^ w^, and thus 



represents a NLO correction. 

The NRQCD calculation of the tt cross section 
sums up the insertions of Coulomb potentials to 
all orders in as- Since the suggested OPE ex- 
pansion of the phase-space corrections necessar- 
ily involves a matching of diagrams with a given 
order in v and as, we have to test also the con- 
vergence of the matching corrections arising from 
the addition of Coulomb potentials. The inser- 
tion of a Coulomb potential yields a factor as/v 
for A ^ cxd; therefore we shall expect a suppres- 
sion factor asTTit/A for every Coulomb exchange 
added to the matching relation of Fig. [D An op- 
timal range of values for the cutoff is to be de- 
termined in order to find the required suppres- 
sion of both the {h? /rrif)^ power-counting break- 
ing terms and the {asmt/A)^ Coulomb phase- 
space corrections. Let us note that the phase- 
space matching of the (e+e~)(e"'"e~) operators at 
2-loops requires the removal of A-dependent non- 
local terms by 1-loop diagrams with phase-space 
subtractions in the production currents, in anal- 
ogy to the process of removal of subdivergences 
in the usual renormalization procedure. 

Results for all the phase-space corrections 
needed for a determination of the threshold tt 
cross section with NNLL accuracy and a com- 
plete numerical analysis shall be presented in a 
forthcoming publication. 
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